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CN ' Abstract 

o: 

Quantization is performed of a Friedmann-Robertson- Walker universe filled with 
a conformally invariant scalar field and a perfect fluid with equation of state p = ap. 
A well-known discrete set of static quantum wormholes is shown to exist for radiation 
(a = 1/3), and a novel continuous set is found for cosmic strings (a = —1/3), the 
. latter states having throat radii of any size. In both cases wave-packet solutions 

to the Wheeler-DeWitt equation are obtained with all the properties of evolving 
quantum wormholes. In the case of a radiation fluid, a detailed analysis of the 
quantum dynamics is made in the context of the Bohm-de Broglie interpretation. 
It is shown that a repulsive quantum force inversely proportional to the cube of the 
scale factor prevents singularities in the quantum domain. For the states considered, 
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there are no particle horizons either. 
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1 Introduction 



Quantum cosmology and speculations about quantum effects in the very early Universe 
began with the work of DeWitt |l|]. Since then, many investigators have been facing the 
difficult task of giving a reasonable description of the Planck era. Friedmann-Robertson- 
Walker (FRW) models have been intensely studied, and it has been found that some 
of them possess wormhole quantum states. Quantum wormholes are solutions to the 
Wheeler-DeWitt equation which remain regular as the three-geometry collapses to zero 
and are exponentially damped for large three-geometries pf. In quantum cosmology 
wormhole states are not restricted to exotic types of matter fields, and in fact they exist 
in models with a simple matter content, such as dust and a conformally coupled scalar 
field 1. 

Recently, in collaboration with J. C. Fabris and F. G. Alvarenga, we quantized a FRW 
model having as matter content a perfect fluid with an arbitrary barotropic equation of 
state jl. Wave packets were constructed and the behaviour of the scale factor discussed 
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according to the many-worlds and the Bohm-de Broglie interpretations. Since fundamen- 
tal fields are believed to play an essential role in the dynamics of the early Universe, in 
the present paper we study the consequences of adding a conformal scalar field to the 
matter content. The treatment relies on Schutz's canonical formalism ||, which describes 
a relativistic fluid interacting with the gravitational field and endows the fluid with dy- 
namical degrees of freedom. The quantum properties of the model are investigated on 
the basis of the associated Wheeler-DeWitt equation. The super-Hamiltonian constraint 
turns out to be linear in one of the momenta, and a time variable can be naturally in- 
troduced, reducing the Wheeler-DeWitt equation to a bona fide Schrodinger equation. 
The presence of the conformal scalar field gives rise to significant changes in the quantum 
dynamics. A well-known discrete set of static wormhole wave functions exists in the case 
of radiation. In the case of cosmic strings a new continuous family of quantum wormhole 
states is found whose throat radii can be arbitrarily small or arbitrarily large. Both for 
cosmic strings and radiation, wave packets are obtained with the properties of evolving 
quantum wormholes. For a radiation fluid we manage to construct the exact propagator, 
which allows us to study the evolution of initial states of the Gaussian type. We show 
that a strong repulsive quantum force, which is inversely proportional to the cube of the 
scale factor, prevents the appearence of singularities. There is no particle horizon either. 

This paper is organized as follows. In Section 2 a Hamiltonian treatment is given 
to a FRW universe filled with a perfect fluid (with an arbitrary barotropic equation 
of state p = ap) and a conformally invariant scalar field. In Section 3 the classical 
equations of motion are solved for some particular values of a. In Section 4 the Wheeler- 
DeWitt equation is constructed. Static quantum wormholes are obtained in Section 5 for 
radiation and cosmic strings. In Section 6 dynamic wormholes are obtained for cosmic 
strings and radiation, in the latter case by means of the exact propagator to the Wheeler- 
DeWitt equation. The rest of the paper deals only with radiation. The Bohm-de Broglie 
interpretation is employed in Section 7 to study the quantum dynamics, and the role of 
the quantum force to prevent singularities is highlighted. Section 8 is dedicated to final 
comments. 



2 FRW Model With Perfect Fluid and Conformal 
Scalar Field 



A homogeneous and isotropic cosmological model is characterized by the Friedmann- 
Robert son- Walker metric 



where cr^ denotes the metric for a 3-space of constant curvature A; = +1,0 or— 1, and the 
lapse function has been conveniently parametrized as N(t)a(t). 

In units such that c = IQ-kG = 1, the pure gravitational action is 



ds 



,2 



N(t) 2 a(tfdt 2 + a(t) 2 a ij dx i dx j 





(2) 
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where K is the trace of the extrinsic curvature Ky of the boundary dM of the space-time 
manifold M. The matter content is a perfect fluid plus a scalar field conformally coupled 
to gravity. The action associated with the sources of gravity is 



S m = J d 4 x^ p-l [ d A x^ (d^dy + - i- / d 3 xVh K<\> 2 . (3) 

JM 2 JM V 6 / 12 J dM 

Schutz's canonical formalism [[| makes use of a representation for the four-velocity of 
the fluid as U v = - (e,„ + (j3, u + 9S, U ) , where eX,/3,0,S are five velocity potentials. 
The specific enthalpy /i is expressed in terms of the velocity potentials by means of the 
normalization condition U V U V = — 1. The potential S is the specific entropy, and in FRW 
models the potentials ( and (3 are zero. 

Compatibility with the homogeneous spacetime metric is guaranteed by taking the 
scalar field and the velocity potentials as functions of t only. Taking p = a p as equation 
of state for the fluid, and performing an ADM reduction described in detail in [§,0, we 
can write the total action St — S g + S m as || 0, |8| 

S T = j dt(p a a + + p,e + p s S - NH) , (4) 

where we have used the constraint pe = and redefined the scalar field by means of 
$ = a<f). The super-Hamiltonian 7i is given by 

U = -4 - Qka 2 + ^ +1 a 1 - 3 " e s + 4 + 6A:$ 2 . (5) 
24 24 v ' 

The canonical transformation [Q 

T = p s e~ s p- e ^ +l \ p T = e s P : +1 , e = e -(a + l)^, p- e = p € (6) 

Pe 

reduces the super-Hamiltonian to 

2 2 

H = -^ + ^-6ka 2 + 6k<S> 2 + a 1 - 3a p T . (7) 

Since e and pi do not appear in the super-Hamiltonian (0), they do not represent a 
dynamical degree of freedom and may be simply dropped from the action, which takes 
the reduced form 



S = J dt{p a d + p^ + p T f - NH) 



3 The Classical Equations of Motion 

The variational principle 5S = leads to the classical equations of motion 
d = -^ Pa , $ = ^$ , f = Na 1 ~ 3a , 

(9) 

p a = N [I2ka - (1 - 3a)a- 3a p T ] , p^ = -l2Nk$ , p T = 
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supplemented by the super-Hamiltonian constraint 



2 2 

_ Pa + P* _ Q ka 2 + Q k($) 2 + fl l-3a = Q . ( 10 ) 
24 24 ^ v 7 

It follows that pt = constant, and in the conformal-time gauge (N = 1) we have 

$ + M> = , T = a 1 " 3 " , a = -fca + ~ 2 a a~ 3a p T . (11) 

The cases a = — 1 and a = have been previously studied in ||. If a = 1/3, 
the equation of motion for the scale factor becomes identical to that satisfied by $ . If 
a 7^ 1/3 , equation ( ]TT| ) for the scale factor corresponds to Newton's equation of motion 
for a particle of unit mass under a force whose potential is V(a) = ka 2 /2 + ia}~ ia , where 
t = —pt/12 is a constant. Making use of "energy" conservation we find 

r da 

t ~ J V£ - ka 2 - laF^ ' (12) 

where £ is an integration constant. The above integral can be expressed in terms of 
elementary functions in the cases of radiation (a = 1/3), stiff matter (a = 1) and cosmic 
strings (a = — 1/3). In other physically interesting cases such as dust (a = 0), domain 
walls (a = —2/3) and vacuum (a = —1) the integral is given by Jacobi elliptic functions. 

For future reference we solve the classical equations of motion in the radiation case 
(a = 1/3). In the conformal-time gauge (N = 1) eqs.@ are solved by (the three lines 
below correspond to k = — 1, 0, 1, respectively) 

o(r) = Ci sinh r + C2 cosh r , $(r) = c 3 sinh r + C4 coshr , 

o(r) = c\T + c 2 , &(t) = c 3 t + , (13) 

a(r) — C\ sin r + c 2 cos r , $(r) = c 3 sin r + c 4 cos r , 

where Ci, c 2 , c 3 , c 4 and are integration constants. For k = 1 both a(r) and $(r) 
are oscillating functions. The cases k = — 1 and k = correspond to an ever-expanding 
Universe. 



4 Quantization 

The Wheeler-DeWitt quantization scheme consists in promoting the canonical mo- 
menta to operators according to 

Pa - -i^ , P. - , PT - -i^ , (14) 

and, with the corresponding super-Hamiltonian operator TC, forming the Wheeler-DeWitt 
equation 

m?(a, T) = 0, (15) 
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where \l/(a, $,T) is called the wave function of the Universe. In our case, the Wheeler- 
DeWitt equation associated with eq.fllOD takes the form of the Schrodinger-like equation 



For the sake of convenience, let us perform the reparametrization 

R - x 



which casts equation ([H]) into the form 



T 



-T 



(17) 



1 a 2 



2dR 2 2 



1 a 2 fc , 
h — Y" 

2 9x 2 2 A 



i? x 1 - 30 dV(R,x,r) 



dr 



Eq.flTSD can be written as id^/dt = H^S? with a self-adjoint Hamiltonian operator if the 
inner product is [ffl 



/oo /*oo 
-oo JO 



(19) 



5 Static Quantum Wormholes 

For the radiation case (a = 1/3) and k = 1 the stationary solutions of the Wheeler- 
DeWitt equation (|18|) take the form ^(R,x,t) = ip n (R)(p m (x) e ~ lET where <p n is the n-th 
eigenfunction of the harmonic oscillator and E = E nm = n — m, with n, m — 0, 1, 2, . . .. 
The wave functions ip nm (R,x) = l Pn{R) l Pm{x) form a well-known set of wormhole quan- 
tum states g g,|nj|nj. 

We now show that static wormholes also exist for cosmic strings and k — 1. With 
a = —1/3 the Wheeler-DeWitt equation flTgp becomes 



i # 



1 



2d X 



2 f ^ X 2 )vI/(i?, X ,r)=^^ 7 ^-vl/(i ? , x , r ). (20) 



I? d_ 
12 «9r 



Writing = e ip we find 



<9 2 </> 

aR 2 



+ kR 2 i) 



<9 2 V 
<9x 2 



12 r 



(21) 



Applying the method of separation of variables with if>(R, x) — Y(x) we § e t 

^- x " + kR2x ^ E T2 = ^- Y " + k ^ = l' 

where A is a separation constant. Thus we have the following ordinary differential equa- 
tions: 



:22i 
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X" + (A-kR 2 + E —)X = , (23) 

6 

Y" + {A - k X 2 )Y = . (24) 
With k = 1 and A = these equations reduce to 

X"- - j^jR 2 X = , (25) 

Y"- x 2 Y = . (26) 

The general solution to these equations which is regular everywhere and exponentially 
damped for both large scale factor and large scalar field is 



X(R) = R 1/2 K 1/4 (<Jl - E/6 R 2 /2) , (27) 

Y( X ) = |x| 1/2 ^i/4(x 2 /2) , (28) 

where E < 6 and K v is the modified Bessel function of the second kind, of order v. 
Thus, the static wormhole wave functions in this case are 

MR,X) = N E (| X | J R) 1 / 2 K 1/4 ( X 2 /2)^i/4( sfl ~ E/6R 2 /2) , (29) 

with E < 6. These wave functions are normalizable in the inner product ©. The 
normalization constant is found to be [Ref . |i~3|] , formulas 6.521(3) and 6.576(4)] 



7T 



N e = Z^\/ 1 - E / Q ■ ( 3 °) 



Thus, in the case of a fluid of cosmic strings there exists a continuous family of wormhole 
quantum states. Because the Hilbert space of states is separable, the wave functions (|29| ) 
cannot be mutually ortogonal, just like the elements of the overcomplete set of coherent 
states of the harmonic oscillator. Indeed, we have 

( , , ] = 4(l-^/6) 3 / 8 (l- W /8 r , n 

We>, We) ^ _ E ,/ 6 y /4 + ^ _ £/ 6 )i/4] [(! _ E'/Qy/ 2 + (1 - E/6)V 2 } ' 1 ' 

The expectation value of the scale factor in the state fl29| ) is 

(R) - T{1/4)2 (32) 

{R) ~ (1-E/6)V4 ' ^ 32 ) 

which may be interpreted as the radius of the throat. Since — oo < E < 6, static worm- 
hole universes with arbitrarily small or arbitrarily large throat radii are allowed. From 
equations (0), ([H]) and \l/ = e~ lEr ip it follows that pt = E. For E < 6 the classical 
equations of motion ( ITT] ) with a = — 1/3 have no static solution for the scale factor. Thus, 
the quantum regime in the state (|29|) is such that the classical evolution of the universe 
is halted. 
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6 Evolving Quantum Wormholes 



In the case of cosmic strings, let us consider E > 6. Setting E = 6(1 + e 2 ) with e > 0, 
a solution to eq.(|25|) is 



X(R) = R l/2 J l/i (eR 2 /2) , (33) 
where J v is the Bessel function of the first kind, of order v. The stationary wave functions 

A(R,X) = (\ X \R) l/2 K 1/4 tx 2 /2)JideR 2 /2) (34) 
do not represent physical states because they are not normalizable in the inner product 
(19|). Normalizable wave packets may be constructed by writing 



C(e)e-^ £ (R,x)de 



e-^(\x\R) l/2 K 1/4 ( X 2 /2) I C(e)e- 6 ^ 2 J 1/4 (eR 2 /2)de 



The choice 



C(e) 



e 5/4 e - 7 ^ 



with 7 a positive constant, leads to [Ref . |1| , formula 6.631(4)] 



*{R,X,t) 



-6ir 



(7 + 6ir) 5 / 4 



\X 



1/2 ir 1/4 (x 2 /2)i?exp 



I?. 1 



16(7 + 6ir) 



(35) 



(36) 



(37) 



is 



(3f 



up to a normalization factor. The expectation value of the scale factor in the state 

W(T) " rwi) I - —) ' 

Note that the mean radius of the Universe grows without bound in spite of the fact that 
the spatial geometry of the Universe is closed (k = 1). Taking 7 sufficiently small, so that 
the wave packet is very well localized near R = at r = 0, the expectation value of the 
scale factor becomes as large as one pleases immediately after r = 0. In other words, it is 
nearly certain that the Universe will have an enormous radius at any positive time. This 



was called "inflation without inflation" in 14, 15 



Let us take up again the radiation case. With a = 1/3 we can put equation ( jIB|) in 
the form 

i S-*{R, x, r) = H#{R, x, r) = (H R - H X )*{R, X , r) (39) 



dr 



where 



H 



1 d 2 k 2 
R ~ 23R 2 + 2 



1 d 2 



k 



Hx ~ 2d X 2 + 2 r 



The Hamiltonian operator H is self-adjoint in the inner product 

/oo roo 
d X dRV(R, X y$(R,x) 
-00 Jo 



(40) 



(41) 
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as long as its domain is restricted to those wave functions such that \I/(0, x) — or 
^'(O, x) = 0, where the prime means partial derivative with respect to R ||15|| . Except 
for this domain restriction, H is the difference between: Hamiltonians of two harmonic 
oscillators (k = 1); Hamiltonians of two free particles (k = 0); or Hamiltonians of two 
inverted oscillators (k = —1). 

The propagator for ( P§| ) is 

G( X , R; x', R'\ r) = Gr(R] R'; t)G*( X ; x'] t), (42) 

where G is the usual harmonic-oscillator propagator. As concerns Gr, some care must 
be exercised since R takes values only on the half-line R > 0. For the sake of simplicity, 
we consider only those wave functions that satisfy the condition ^/'(O, Xi T ) = 0- F° r this 
boundary condition, the propagator is given by |TB| 



G R (R, R'; r) = G(R; R'; r) + G(R; -R'; r) 



where 

G(x, x'; t) 



Vk 



1/2 



Am sin Vk 



exp< 



iy/k 



2 sin y/k t 



[x 2 + x' 2 ) cos Vk t — 2xx' 



(43) 



(44) 



is the propagator for a one-dimensional harmonic oscillator with mass m = 1 and fre- 
quency w = Vk. 

The quantum dynamics of the models will be studied by following the time evolution 
of the initial Gaussian wave function 



71 



where a and (3 are positive constants. Making use of 

/oo /*oo 
dx! dR'G( x ,R;x',R';r)^(x\R',0) 
-oo JO 

the wave packet takes the form 



(45) 



(46) 



'4v^5 



7T 



\ cos 2 (y / A;r) 2a tan(Vkr) — iVk 2(3 tan(Vkr) + iVk 



x exp 



2 tan(VkT 
i\fk 



iVk 



2tan(v^r) 



1 + 



cos 2 (\/A;t) 2/3tan(A/fcr) + 



cos 2 (-\/A;t) |^2crtan(v / A;T") — 



I? 2 



(47) 
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7 The Bohm-de Broglie Interpretation 



The Bohm-De Broglie interpretation ]T(J is an alternative interpretation of quantum 
mechanics which allows the treatment of a unique system such as the Universe, which 
cannot be repeatedly prepared in the same state as required by the standard statistical 
interpretation. The Bohm-De Broglie interpretation starts by writing the wave function 
of the Universe in the form 

= 0exp(iS), (48) 

where and S are real functions. Then, inserting (|48"D into the Wheeler-DeWitt equation 
(|l6|) with a = 1/3, there results 



OS, 



1 
2 



dS\ 2 _ fdS\ 



dS 



+V+Q=0 



fdQ\ (dS 
VUR VUR 



+ 



d 2 s 



d x 2 



d 2 s 
W 2 



(49) 



where V and Q are the classical and quantum potentials, respectively: 



V 



-k(R 2 

9 V 



X' 



Q = 

V 20 



d 2 e o 2 e 



dR z dx 1 



(50) 



Let us first examine the case of vanishing spatial curvature (k = 0). 



7.1 The Flat Case 



The wave packet is obtained by taking the limit as k — > of the wave function (f47|) , and 
it takes the form 



where 



Mx,t)= J 



^(3 



7Y 



*i(i2,r) Mx,r) 



(51) 



2ar 



cxp 



2ot 



2t 2ot 



R' 



(52) 



1 



2(3r + i 



exp 



i 2pr s 
27 2(3r + i 



X 



The wave function (|51|) may be written in the form (|48"D with the functions and S 
given by 

2a 2 r _ Q 2/3 2 r 2 

X +/o(r), (53) 



5 1 + 4<x 2 r 2 
© = 00 CO exp < 



1 + A(3 2 t 2 
aR 2 (3 X 2 



{ 1 + Aa 2 r 2 l + 4/3 2 r 2 )' ^ 

where / an d go are functions of r alone that play no role in the subsequent discussion. 
From this we can find the Bohmian trajectories for the scale factor a(r) and the scalar 
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field 3>(t), recalling that a, R and x are related by flT7|), and also that in the r-gauge 
one must take N = — 1 in Hamilton's equations (^). The Bohmian trajectories are the 
solutions of the differential equations 

Pn = M , Px = ^- (55) 

Taking into account the previous remark on the lapse function, the classical equations of 
motion and the relations (0), we get 

R ^ = YV^ R{T) - (56) 

An immediate integration yields 



R(r) = A Vl + 4a 2 T 2 , (57) 

where Aq is a positive constant of integration. There is no singularity since the Bohmian 
trajectories R(r) never reach R = 0. 

In a similar fashion, we find for the scalar field 

A8 2 t 

*M = i^g5^x(T), (58) 

whose solution is 



X (r) = B ^l+A/3 2 r 2 , (59) 

where B is a real constant of integration. The probability distribution of x derived from 
(|5j]) is an even function of x- Thus, averaging over Bq we obtain the expectation value 
(x) = 0, as it should. 

The quantum potential Q is found by inserting (Q) into (|50|): 

n(R g 2 ° 2R2 2 ? 2 X 2 (an, 

WX,TJ l + W 1+4/5^ (l + 4a2 r 2) 2 + (1 + 4/3V2) 2 ' lDUJ 

The quantum force associated with the time evolution of the scale factor is 

tR ~ dR ~ (I + Aa 2 r 2 ) 2 ■ ^ 

With the help of the trajectory fl5"7D , the quantum force F R can be expressed as a function 
of R only: 

Fr\R) = §- 3 > (62) 

where C is a positive constant. This force is always repulsive away from R — , is 
strongest when -R(t) is minimum and is responsible for the avoidance of the singularity 
in the quantum domain. 
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7.2 The Case of Positive Spatial Curvature 

Now the wave function (E7|) takes the form 



*tox,T) = tP^eexp(iS) , 



7T 



(63) 



with 5 and given by 

© = gi(r)exp 



oR 2 



Px 2 



cos 2 t + 4<t 2 sin r cos 2 r + 4/5 2 sin r 



(64) 



and 



^ = h(r) 



R 2 



2tanr 



1 - 



1 



cos 2 r + 4cr 2 sin 2 r 



r 



2tanr 



1 



cos 2 t + 4/3 2 sin r 



(65) 



The Bohmian trajectories are determined by 



R(t) 



X{r) 



tanr 



tanr 



1 



cos 2 t + 4<t 2 sin r 



1 - 



cos 2 r + A(3 2 sin r 



i2(r) 



X(r) 



(66) 



(67) 



Integration yields 



R(t) = A + \J cos 2 r + 4cr 2 sin 2 r , x( r ) = B + J cos 2 r + 4/? 2 sin 2 r , 



(6* 



where A + and -B+ are constants of integration, with A + > 0. Both R and \ oscillate 
eternally between minimum and maximum values. 
The quantum potential is 



Q(R, X ,r) = h(r) 



2a 2 R 2 



+ 



(cos 2 r + 4cr 2 sin 2 (cos 2 r + 4/5 2 sin 2 



(69) 



from which it is easily seen that the quantum force is again of the form 

Only the interval [0, 7r] is physically acceptable for r at the classical level, correspond- 
ing to the initial and final singularities. In the quantum realm, however, r may take values 
over the whole real line, there being no singularity in virtue of the repulsive character of 
the quantum force. Differently from the flat case, in which only the quantum potential 
exists, in the present case [k = 1) the classical an quantum potentials compete. However, 
as equation (|62"D shows, the repulsive quantum force is much stronger that the classical 
force near the singularity. 
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7.3 The Case of Negative Spatial Curvature 



The results for k = — 1 follow from those for k = 1 by simply replacing the trigono- 
metric functions by the corresponding hyperbolic functions. The Bohmian trajectories 
satisfy 



R(r) 



1 



tanhr 



1 - 



cosh r + 4<t 2 sinh r 



R(r) 



(70) 



X{r) 



which are solved by 



1 



tanhr 



1 



cosh r + 4/3 2 sinh r 



X{r) 



R{t) = A- V cosh 2 t + 4(x 2 sinh 2 r , x( T ) = B - V cos h 2 r + 4 /3 2 sinh 2 r . 
with A_ > 0. Once again the quantum force as a function of R is given by (|62| 



(71) 



(72) 



8 Conclusion and Final Remarks 

In this work we quantized Friedmann-Robertson- Walker cosmological models having 
for matter content a perfect fluid and a conformal scalar field. We found static wormhole 
wave functions in the cases of radiation and cosmic strings. In the radiation case the 
discrete set of wormhole quantum states was previously known in the literature. In the 
case of cosmic strings we found a new continuous family of static quantum wormholes 
with arbitrary throat radii. The members of this continuous set of finite-norm wave 
functions are not orthogonal to each other because the Hilbert state of states is separable. 
In this aspect they resemble the overcomplete set of coherent states of the harmonic 
oscillator. For cosmic strings an evolving-wormhole wave packet was constructed leading 
to a dynamically open universe, although the spatial geometry is closed. If the initial 
state of the Universe is sufficiently localized near the classical singularity, the scale factor 
is enormously large immediately after the initial instant, a phenomenon described in the 
literature as "inflation without inflation". In the radiation case, with the help of the 
exact propagator we constructed wave packets that are examples of dynamic quantum 
wormholes. The behaviour of the scale factor was studied according to the Bohm-de 
Broglie interpretation of quantum mechanics. In all three cases (k = 0, ±1) the singularity 
is prevented by a quantum force inversely proportional to the cube of the scale factor. 
Also, in all three cases, the integral jl^ dr' / R(r') is divergent, so that for the states 
considered the models have no particle horizon. 
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